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There are several ways to explain the dark matter relic density other than by the ordinary freeze-out scenario. 
For example, the freeze-in mechanism may constitute an alternative for generating the correct relic density for 
dark matter candidates whose predicted freeze-out abundance is too low due to a large total annihilation cross 
section. Here we show that although such a mechanism could explain why a dark matter candidate has the 
correct relic density, some candidates may still be ruled out because they would lead to a large gamma ray flux 
in dwarf spheroidal galaxies or a large elastic scattering rate in direct detection experiments. To investigate this 
scenario we examine neutralino dark matter in the MSSM. However our conclusions can be generalised to other 
types of annihilating DM candidates with a low relic density in the freeze-out scenario but which have their relic 
densities generated by some other mechanism. 



I. INTRODUCTION 

For several decades, the vast majority of dark matter (DM) 
models have assumed that DM exists in the form of annihilat- 
ing particles whose relic density is determined by the freeze- 
out mechanism | ljt3J. This mechanism intimately links the 
total DM annihilation (or co- annihilation) cross section to the 
DM relic abundance and leads to the prediction that the total 
annihilation cross section must be of order the weak scale. 

In some models, nevertheless, the total annihilation cross 
section at chemical decoupling can be so large (notably if it 
involves annihilations through resonant channels) that the pre- 
dicted abundance for the DM candidate is well below the ob- 
served value (41, ruling it out as the dominant contribution to 
the DM relic abundance. However, new mechanisms, such as 
the decay of a metastable particle or the freeze-in scenario Q, 
have been proposed as viable alternatives to restore the relic 
density of such candidates to the required value. 

The freeze-in mechanism can be relevant for models where 
a feeble coupling exists between the thermal bath in the early 
Universe, consisting of visible sector particles (e.g. standard 
model or MSSM particles etc), and a thermally decoupled par- 
ticle species, hereafter referred to as Feebly Interacting Mas- 
sive Particles (FIMPs). Such a coupling is too small to main- 
tain the FIMPs in thermal equilibrium but it can nevertheless 
lead to their production via the freeze-in mechanism |5|. In 
essence, through some interaction or decay of particles from 
the visible sector, energy density leaks out of the thermal bath 
in the form of FIMPs as once they are produced it is unlikely, 
due to a small number density and the suppressed interactions 
with the visible sector, that they can reproduce visible sector 
states, and so their number density accumulates until the pro- 
cess that produces them drops out of equilibrium |5 1. 

The freeze-in mechanism can lead to a number of possibili- 
ties for DM. For example, if the FIMP is stable (via some sta- 
bilising symmetry) it can itself play the role of DM. Its abun- 
dance is set directly by the freeze-in mechanism and depends 
on the size of the interaction or decay rate that leads to the 
production of FIMP states. 



A second possibility, which we pursue as a working sce- 
nario in this paper, is that these FIMPs are unstable and decay 
into other particles that will constitute the DM. If the FIMP 
transforms under the same stabilising symmetry that keeps the 
DM particle stable and has a larger mass, then decays of the 
FIMP states will generically produce an abundance of DM 
particles. Due to their small coupling the FIMPs will have a 
sufficiently long lifetime such that they will decay after the 
DM freezes-out and regenerates the DM abundance. In sce- 
narios where the freeze-out abundance of DM is too low, the 
decay of the frozen-in FIMPs can reproduce the correct abun- 
danc4B 

The question that we wish to address in this paper is 
whether such a mechanism (or related mechanisms) can both 
explain the observed DM abundance in models where the DM 
annihilation (or co- annihilation |6, 7 |) rate at freeze-out is too 
large, and be simultaneously compatible with the latest results 
from direct and indirect detection experiments. In particular, 
for direct detection constraints, we apply limits on the DM- 
nucleon spin-independent elastic cross section as derived from 
the XENON 100 experiment lO. The indirect detection limits 
that we apply come from the latest observations of the dwarf 
spheroidal galaxies (dSphs) by FERMI-LAT |9 |, which place 
an upper limit on the gamma flux emerging from DM annihi- 
lations. 

We consider as our DM candidate the neutralino of the 
MSSM and assume that one can add an extra term or terms 
to the MSSM Lagrangian in order to implement the freeze-in 
scenaric]^ The use of Supersymmetry enables us to explore 
very different types of configurations in terms of resonances 
and co-annihilations. However, similar conclusions will also 
hold for other types of DM candidates where the freeze-out 



An appropriate choice for the size of the feeble coupling is needed to get the 
correct abundance. We treat this coupling as a free parameter and assume 
that it can be chosen correctly such that the correct final DM abundance is 
generated. 
^ See ^JO|JJJ for examples. 
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Scan A 
2 GeV < Ml < 120 GeV 



ScanB 
90 GeV < Ml < 2000 GeV 



90 GeV < M2 < 2000 GeV 



200 GeV < Ms < 6000 GeV 



2 GeV <]u< 2000 GeV 



0.1<tanp<75 



-4000 GeV <At< 4000 GeV 



100 GeV < m^o < 1500 GeV 100 GeV < m^o < 4000 GeV 



100 GeV <mr < 4000 GeV 



100 GeV < < 4000 GeV 



100 GeV < m^^ 2 < 4000 GeV 



100 GeV < < 4000 GeV 



TABLE L Allowed ranges of the parameters. 



relic density is too low due to a large annihilation cross sec- 
tion. 

In Section we describe our method for investigating the 
parameter space of the MSSM with the aim of determining 
the regions that lead to an annihilation cross section greater 
than the common value of Gv ~ 3 x 10~^^cm^s~^, which is 
needed in the standard freeze-out scenario. This is indeed 
where freeze-in (and similar mechanisms) can be important 
as regions of parameter space that were previously unable 
to explain why the DM candidate would make up all of the 
WMAP observed value are now potentially viable. In Sec- 
we investigate the phenomenology of these "under 



tion 



III 



abundant" configurations in the light of DM experiments and 
conclude in Section Hvl 



II. PARAMETER SPACE 

We consider the MSSM and allow for 11 parameters to 
vary, namely the gaugino masses, Mi, M2, M3, the Higgs- 
Higgsino mass parameter, ju, the ratio of the Higgs vacuum 
expectation values, tanP, the stop trilinear coupling. At, (all 
other trilinear couplings are set to zero), the mass of the CP- 
odd Higgs, m^o, and finally the parameters m^^2' ^gs 
m^~^^, which represent the squark masses for the first two gen- 
erations, the third generation squark masses, and all genera- 
tions of the "left" and "right" sleptons respectively. 

The choice of m^^ = rriq^ ^ ruq^ is particularly relevant 
since stops can be lighter than the first two generations and 
can be relevant for enhancing neutralino annihilations (cf e.g. 
lfT2l ). Separating the "left" and "right" slepton masses also al- 
lows for a light slepton (mostly in the case of "left" sleptons) 
that can play a significant role in neutralino co-annihilations 
(131 [T4I. The values of all parameters are defined at the elec- 
troweak scale. 

We perform two Markov Chain Monte Carlo (MCMC) 
scans labelled Scan A and Scan B. Scan A is dedicated to 
low neutralino masses (below 100 GeV) while Scan B is ded- 
icated to heavier candidates (above 100 GeV). This choice of 
two separate scans above and below 100 GeV is purely arbi- 
trary but it turns out to be a useful division. The reason being 



Constraint 


Value 


Tolerance 




< 0.1123 |15 | 


none 




25.5 X 10-^0 


Stat: 6.3 X 10-10 
sys: 4.9 x 10-^^ 


Ap 


< 0.002 


0.0001 


B¥{b -> si) 


3.55x10-4 1161 


th: 0.24 X 10-4 




< 4.5 X 10-9 (ni 


4.5 X 10-11 


R(B -> TVi:) 


1.36 [H 


0.23 


r(z^x?x?) 


< 1.7 MeV 


0.3 MeV 




<0.1pb |18| 


O.OOlpb 



TABLE II. Constraints used to calculate likelihoods, from Ref. [191 
unless stated. Here Ilpo^^ is the relic abundance of neutralino DM 
from freeze-out, Ap is the contribution to the electro-weak precision 
variable p, R(5 xv^) is the ratio of the MSSM to SM branching 
fraction of 5+ T+Vx- 



that the neutralino candidates found in each scan represent dif- 
ferent freeze-out scenarios and are most sensitive to different 
experimental searches. Scan A features s-channel resonant 
effects while Scan B shows a greater number of t-channel ex- 
change and co-annihilation processes. In these scans we are 
looking for points where the freeze-out relic density, Q^mh^, 
is lower than the mean value ^^wmap^^ 0.1123 (obtained 
by combining WMAP data with BAO and measurements 
1 15 1) and which satisfy a number of constraints from particle 
physics experiments. We take the mean value for ^^wmap^^ as 
an absolute limit with no uncertainty rather than the WMAP 
maximal upper bound as we are interested in only those points 
that would require regeneration in order to fit the WMAP data. 
Since the most interesting region of the parameter space in 
terms of regeneration is far away from the limit, the uncer- 
tainty in the mean value can be neglected without affecting 
the results. Our choices for the allowed ranges of the MSSM 
parameters for the two scans are listed in Table |l] and the con- 
straints used to calculate likelihoods can be found in Table [III 

All the physical quantities in this analysis are computed us- 
ing the micrOMEGAs code 1 20 1 except for the SUSY particle 
spectrum and decay rates of the Higgs particles which were 
calculated using SoftSusy 1211 and SUSYHIT |22| respec- 
tively. LEP limits on the sparticle masses are applied automat- 
ically by micrOMEGAs (see |20| for details). In addition, a 
lower limit on the squark masses is set at 100 GeV. The Higgs 
masses are restricted to the allowed range by the HiggsBounds 
programme [2124] . SoftSusy, SUSYHIT, micrOMEGAs and 
HiggsBounds were interfaced via the SUSY Les Houches Ac- 
cord 1 251. We do not apply limits on the squark and gluino 
masses coming from the latest CMS and ATLAS data. The 
effect of these limits could be considered by simulating events 
for a converged sub set of the Markov chains as described in 
1261 , however, we consider this beyond the scope of this par- 
ticular work. Our focus here is to examine the possibility of 
regenerating the DM density in under- abundant DM scenarios 
in the light of DM experiments. 

To explore the parameter space we generate a random walk 
using the Metropolis-Hastings algorithm. For these scans a an 
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Standard model parameter 


Mean value 


Experimental uncertainty 


rut 


172.9 GeV 


1.5 




4.19 GeV 


+0.18 
-0.06 




0.1184 


0.0007 




127.916 


0.015 



TABLE III. Constraints used to calculate likelihoods for standard 
model parameters, from Ref. 1 19 |. 



initial point in parameter space is randomly chosen. Follow- 
ing this steps in the random walk are taken along randomly se- 
lected directions in the parameters space and an initial "burn- 
in" phase is used to adjust the magnitude of the proposed step 
size for each direction to optimise the exploration of the pa- 
rameter space, this is periodically adjusted during "burn-in" to 
ensure that the parameter space is covered as fully as possible. 
The directions in which steps are taken were generated from 
the eigenvectors of the covariance matrix found in preliminary 
scans. The "burn-in" phase also ensures that the chain has al- 
ready converged towards a high likelihood before points are 
recorded. The total likelihood function is formed by the prod- 
uct of partial likelihoods for each observable in Table [ll] As in 
Ref. |27 1 we use a Gaussian distribution for observables with 
a preferred value 

F2{x,jLi,a) = e 2.2 ^ (1) 

where ju is the preferred value of the observable and a is the 
tolerance. For observables with only an upper or lower limit a 
distribution of the form, 

F3 (x, A/, a) = , (2) 

is used. Here a is positive for lower bounds and negative for 
upper bounds. For the relic abundance, the masses of the spar- 
ticles and the Higgs masses, the partial likelihood is either one 
or zero as no uncertainties are included. 

Uncertainties in standard model parameters were in- 
cluded in the form of nuisance parameters which are then 
marginalised as part of the random walk. The mean values 
and uncertainties of the nuisance parameters are shown in Ta- 
bleini 

A. Scan A: results for scenarios with m^o < 100 GeV 

Previous supersymmetric parameter scans either looked for 
scenarios with the correct relic density (e.g. |28-35|) or re- 
laxed the constraint on the relic density, allowing for very 
small ^Ipo^^, and did not assume the presence of regenera- 
tion mechanism [ 27, 36 1. In this paper we will both relax the 
lower bound on the relic density and assume that the freeze- 
in mechanism can regenerate the relic density to the observed 
value. 

In FIG. [T] we show the relic density versus DM mass for 
candidates found by the MCMC. In most scenarios more than 




102 

m^o [GeV] 

FIG. 1. Plot of Q^Y^oh^ against m^o. The colour coding represents 
the process with the largest contribution to the neutralino annihila- 
tion rate, which determines the freeze-out relic abundance. Green 
points correspond to resonant annihilation via Z, red points to res- 
onant annihilation via the light Higgs boson (h^), orange points to 
resonant annihilation via the pseudo-scalar Higgs (A^), blue points to 
stau CO- annihilation or annihilation via stau exchange, violet points 
to chargino co-annihilations or chargino exchange, black points to 
squark co-annihilation (all squark flavours). 



one process will contribute to the freeze-out relic abundance 
but in FIG.[T]the largest single contribution to the annihilation 
rate, which in the majority of scenarios dominates the others, 
is indicated. In all of the following plots the points found by 
the random walk are plotted as semi-transparent dots, faint re- 
gions therefore correspond to a low density of points while 
regions of strong colour correspond to denser regions. As 
expected there are two visible resonance regions |6|, corre- 
sponding to Z gauge boson and light CP-even Higgs (h^) s- 
channel resonances. In addition there are the usual points 
corresponding to heavier neutralinos that can annihilate via 
s-channel exchange of the CP-odd Higgs (A^) ISTIl , as is well 
known from traditional freeze-out scenarios. These points ap- 
pear as a smeared out region due to the large variation in the 
value of m^o . 

In addition to the s-channel processes the well known t- 
channel exchange and co-annihilations processes involving 
charginos, staus and squarks are also found by the MCMC. 
It is likely that the majority of the points corresponding to 
squark exchange and co-annihilation will be excluded by the 
LHC or Tevatron. However, we still include these points as 
our focus here is to examine the effect of regeneration and the 
resulting DM detection constraints on the possible regions of 
the parameter space. 

The composition of the neutralino LSP in terms of the weak 
eigenstates, the Bino, Higgsinos and Wino differs slightly for 
the various regions displayed in FIG. [T] 

For the Z and resonance regions the neutralino is mostly 
Bino with a small Higgsino component. As is well known, 
(see for example |[37l|38|), the size of the Higgsino compo- 
nent will play a central role in determining the cross section 
for DM annihilations via s-channel Z and h^. This Higgsino 
component will also lead to the dominant contributions to the 
spin-independent elastic scattering cross section in direct de- 
tection experiments, where the main process is the t-channel 
exchange of a Higgs. This connection is important for what 
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m^o [GeV] 



FIG. 2. Plot of Ilpo^^ against m^o . Colour coded for the process with 
the largest contribution to the total neutralino annihilation rate, which 
determines the freeze-out relic abundance. Red points correspond to 
chargino co-annihilation, green points to annihilation via chargino 
t-channel exchange, blue points to annihilation via s-channel Higgs 
(roughly speaking the blue points above the green band correspond 
to annihilation via an s-channel into tt and bb, the few below are 
s-channel annihilation via A^), yellow points correspond to a either 
squark co- annihilation or gluino-gluino annihilations (the latter in 
the case where the gluino is approximately mass degenerate with the 
neutralino DM and its freeze-out sets the neutralino relic abundance). 



relic abundance scales as the inverse of the thermally averaged 
cross section, which in turn scales approximately as the in- 
verse of the neutralino mass squared. As a result the minimum 
relic abundance will increase quadratically with the mass of 
the neutralino. Co-annihilation with light stops is expected 
to add a few more points (below the "quadratic" limit) when 
there is a large fine-tuning between the neutralino and the stop 
mass. However, the stop and neutralino self-annihilation cross 
sections both decrease with the mass of these particles and an 
increase in the fine tuning becomes less and less effective in 
compensating for the lack of efficiency of the co-annihilation 
process when the neutralino mass increases. Besides, these 
points become more difficult to find by the MCMC as they 
require smaller variance (i.e. more dedicated searches). 

The compositions of the higher mass neutralinos is more 
varied than the lower mass states. For example, in points 
whose freeze-out annihilation rate is dominated by chargino 
co-annihilation and t-channel chargino exchange the neu- 
tralino DM can be mostly Wino. For points whose freeze-out 
annihilation is dominated by s-channel Higgs processes, the 
Higgsino component of these neutralinos can be much larger 
(even dominating the composition) than that for neutralino 
DM with masses below 100 GeV. 



follows in the later sections. 

In the cases where t-channel exchange and co-annihilation 
processes, involving light SUSY squarks and sleptons, domi- 
nate the freeze-out dynamics, the neutralino can have a much 
smaller Higgsino component. This is because, in contrast to 
the s-channel annihilation processes, the t-channel annihila- 
tion and co-annihilation diagrams can occur for pure Bino 
neutralinos. 



B. Scan B: results for scenarios with m^o > 100 GeV 



In the case of neutralinos heavier than 100 GeV, one does 
not expect any resonance structure in the (m^o,^^FO^^) plane 
since there are no fixed mass neutral particles (such as the 
light CP-even Higg^or Z boson) that can be produced in an s- 
channel resonance. Instead resonant annihilation through 
will appear over a range of different neutralino masses. Non- 
resonant annihilation via the and Z bosons can still produce 
a large enough cross section to reduce the relic abundance 
for masses above 200 GeV. Chargino or squark t-channel ex- 
change and co-annihilations also lead to an enhanced cross 
section but this does not appear as a fixed mass resonance. As 
a result, we find a smooth homogeneous distribution of points 
in the (m^o,^lFO^^) plane, as shown in FIG.|2j 

The most visible trend in FIG.|2]is that the minimum relic 
abundance found by the MCMC increases quadratically as a 
function of mass. This dependence of the relic abundance on 
the mass of the neutralino DM arises due to the fact that the 



^ Although the /z^ mass is not fixed, it is restricted to a narrow range in the 
MSSM. 



III. DM REGENERATION IN THE LIGHT OF FERMI-LAT 
AND XENONIOO LIMITS 

To examine the impact of a possible regeneration mecha- 
nism we apply limits arising in direct and indirect detection 
experiments to the points found by the MCMC. We do so in 
two cases. The first where there is no regeneration and the 
DM density is set by the value determined by freeze-out. The 
second where regeneration of the DM density has taken place 
after freeze-out and has been regenerated to the WMAP ob- 
served value. The limits for direct and indirect detection are 
applied as 95% confidence level exclusions to the points found 
by the MCMC after the scans have completed rather than in- 
cluding these limits in the likelihood calculations. This allows 
the two scenarios to be compared directly using the same set 
of points. 

We look at the effect of regeneration in the planes 
(asi,m^o), (GshO^Foh^), (<^>pp,m^o) and (<J>pp,asi), where Gsi 

is the spin-independent elastic scattering rate, ^Ipo^^ is the 
relic abundance generated by freeze-out only and <I>pp, which 
encodes the "particle physics input" to the total flux of gamma 
rays from annihilating DM in the dSphs. The quantity <I>pp is 
defined as 



^pp 



(GV) 



-^^0 



dN 
dE 



(3) 



where (av) is the thermally averaged cross section for DM 
annihilation, Eq is the minimum threshold energy considered, 
Emax is the maximum photon energy the limit is sensitive to 
and ^ is the gamma ray spectrum averaged over all of the 
different annihilation channels. Neglecting propagation the 
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m-o [GeV] m-o [GeV] 



FIG. 3. Spin-independent cross section versus neutralino mass. The right panel shows the case with the regeneration of the DM reUc density 
to the correct value, the left panel shows the case without. The limit from XENONIOO 1 8 1 as a limits on the spin-independent cross section of 
a DM species with the WMAP observed relic density is shown as a blue dashed line. The yellow points are excluded by XENONIOO, the red 
points are excluded by indirect detection, grey points are excluded by both and green points survive all constraints applied. 



expected flux of gamma rays from a given source reads as 

= <l>pp X /, (4) 

where / is the DM density integrated along the line of sight 
and over the solid angle and sensitivity of the observation. 

An upper limit on the flux and a particular choice of / 
then set an upper limit on <J>pp. In general the upper limit 
on the flux depends on assumptions about the spectral shape 
of the gamma ray source. Choosing the hardest power-law 
model from 1 9 1 gives an upper bound on the photon flux which 
can be divided by / to give a conservative upper bound of 
<J>PP < 7.5 X 10~^^cm^s~^GeV~^ from observations of the 
Draco dSph by FERMI-LAT (91. However, using a com- 
bined analysis of several dSphs places a stronger limit of 
^PP < 5.0 X 10-^0cm^s-iGeV"^L39J. In this case there is 
no single limit on the gamma ray flux and corresponding / 
value, instead the limit on <J>pp is found by Neyman construc- 
tion [40 , 41 1 where each dSph is weighted by its / value. We 
use this combined limit in what follows. 

For each point found by the MCMC the gamma ray spec- 
trum ^ is calculated using micrOMEGAs and integrated 
from 1 GeV to 100 GeV in order to obtain <l>pp. 

In addition to applying constraints from indirect detection, 
we also apply constraints coming from direct detection ex- 
periments. In particular we apply the limits on the spin- 
independent elastic cross section coming from XENONIOO 
ISl . The spin-independent cross section for each point is cal- 
culated automatically in micOMEGAs and we refer the reader 
to 1 20 1 for details. One important point we do note here is that 
we use the default values for the scalar form factors of the pro- 
ton and neutron as set in micOMEGAs In particular we 
use the default value for the strange quark scalar form factors 
as given in ll2Ql as fs'^ = 0.2594. It is well known that this is 
a source of a large uncertainty in direct detection rates, see for 
example 1 32 , 42 , 43] and can lead to a significant change in the 
predicted cross sections. Astrophysical uncertainties can also 
have an impact on the limits applied, see [44, but again we 
do not allow for these uncertainties. 

Finally, for m^o < 50 GeV, the uncertainties on the ex- 
clusion curve, due to the lack of physical knowledge on 



the energy behaviour of the relative scintillation efficiency, 
are important. The latter do not appear in | 8| because the 
XENONIOO collaboration assumed that the uncertainties on 
the relative scintillation efficiency can be well modelled by a 
Gaussian likelihood centred on the Xgff mean value. It was 
not realised that maximising the global likelihood gives more 
weight to the mean (but not necessarily the physical) value of 
Xeff and does not allow the real (physical) uncertainties on Xeff 
to be taken properly into account |46|. Here we continue to 
use the exclusion curve obtained in [8 1 as a guideline to un- 
derstand the effect of regeneration but a more detailed study 
would require the implementation of all these sources of un- 
certainties in the derivation of the direct detection exclusion 
limit. 

In the following subsections we present a series of dou- 
ble panel figures. The plots corresponding to no regeneration 
(freeze-out only contributions to the DM relic density) are dis- 
played in the left panels. The same points are plotted in the 
right panels but now with the DM density regenerated to the 
WMAP observed value. Note that these scenarios are strictly 
identical in the pairs of plots apart from the DM densities used 
to calculate the limits. It should be noted that in the calcula- 
tions for the indirect detection rates, micrOMEGAS 1201 uses 
by default the value of the DM density determined by WMAP 
1 4 1 not the value predicted by freeze-out, which in the majority 
of our cases will be below the WMAP value. In order to cal- 
culate the gamma ray flux for the under- abundant scenarios, 
the square of the scaling factor, r|, needs to be applied, where 
r| = ^^Fo/^WMAP. Similarly, for the under- abundant scenar- 
ios, the limits on the elastic scattering cross section from di- 
rect detection need to be scaled by r\. 



A. Regeneration in scenarios with m^o < 100 GeV 

In FIG. [3] plots of the spin-independent elastic scattering 
cross section against neutralino mass are shown for Scan A. 
In these plots the red points are excluded by the constraints 
from the FERMI-LAT gamma ray limits from dSphs, yellow 
points are ruled out by XENONIOO direct detection searches. 
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FIG. 4. Spin-independent cross section versus the neutralino freeze-out relic density for m^o < 100 GeV. The right panel shows the case with 
the regeneration of the DM relic density to the correct value, the left panel shows the case without. Colour coding is the same as in FIG. [3] 
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FIG. 5. <I>pp vs the spin-independent cross section for m^o < 100 GeV. The right panel shows the case with the regeneration of the DM relic 
density to the correct value, the left panel shows the case without. The limit on <I>pp shown as a blue dashed line is from the combined analysis 
of FERMI-LAT observations of dSphs |39| . Colour coding is the same as in FIG. [3] 



grey points are ruled out by both and the green points are those 
that survive the constraints applied. 

In the left panel of FIG. [3] no regeneration of the DM den- 
sity is assumed; hence r\ can be small. The result is that for 
points with a low freeze-out relic abundance, like those in the 
Z and resonance regions, the elastic scattering cross sec- 
tion can be large, i.e. above the XENON 100 limit as evalu- 
ated for a DM species with the WMAP observed density, and 
still predict a sufficiently low event rate in a direct detection 
experiment to evade the exclusion limits. 

Also visible is a region around and just below m^o ~ 
100 GeV. Comparing with FIG.[T] this region corresponds to 
the scenarios in which chargino co-annihilations and t-channel 
exchange diagrams dominate during freeze-out. 

If we now assume that the DM density is regenerated af- 
ter freeze-out to the observed value, all points above the 
XENON 100 limit are now ruled out, as shown in the right 
panel of FIG. [3] There are a number of points that are still 
allowed, in particular those that appear in the Z and reso- 
nance regions. The reason for this is that if the neutralino DM 
can annihilate via an on-shell s-channel resonance, the size of 
the couplings needed to give a large enough annihilation cross 
section at freeze-out to reduce the DM relic abundance below 
the WMAP measured value, can be smaller. 

The size of the couplings between the neutralino and both 
the Z and is determined by the size of the Higgsino compo- 



nent in the neutralino, which in turn determines the size of the 
spin-independent elastic scattering cross section. This reduc- 
tion in the couplings will therefore allow some of the points in 
the resonance regions to avoid the direct detection limit, pro- 
vided they correspond to points with close to on-shell freeze- 
out annihilations. Despite this, a significant number of points 
are ruled out by direct detection. 

In FIG.|4]we present plots of the distribution of points found 
by the MCMC in the (asi,^^FO^^) plane. Once again, in the 
left panel of FIG. |4] the DM relic density is kept at the value 
predicted by freeze-out and in the right the DM density is as- 
sumed to have been regenerated to the observed value but is 
plotted as a function of the relic density generated by freeze- 
out for each point. 

Different regions of the plots in FIG. |4] can be identified 
and explained in terms of the connection between the an- 
nihilation cross section in the early universe and the spin- 
independent elastic scattering cross section. There are two 
main regions of points corresponding to different types of 
process that dominate the DM annihilation cross section at 
freeze-out, they are, DM annihilation via s-channel Z or 
and DM co- annihilation with another SUSY particle (usually 
the chargino). With reference to the left panel of FIG. |4j the 
points corresponding to s-channel processes are roughly con- 
tained within the green diagonal band and the yellow points 
above. The co-annihilation points are those below the green 



diagonal band. 

Moving from small to large freeze-out abundances (left to 
right in both panels of FIG.|4]but remaining at a constant spin- 
independent scattering rate, corresponds to moving off-shell 
for the s-channel annihilation rate at freeze-out. That is, the 
mass of the neutralino DM is moving away from either mz/2 
or This reduces the overall annihilation rate and there- 

fore increases the freeze-out relic abundance. 

Moving down the plots in FIG.|4]we move to smaller spin- 
independent elastic scattering cross sections with the size of 
the Higgsino component in the mostly Bino neutralino de- 
creasing, which results in smaller couplings to . The DM s- 
channel annihilation cross section at freeze-out also decreases 
with the decreasing couplings and that effect translates into 
the diagonal slope that can be seen in both plots of FIG. |4] 
The maximum size of the annihilation cross section at freeze- 
out, when the s-channel resonance is on shell, decreases with 
decreasing Higgsino component. Consequently the smallest 
possible value of the freeze-out relic abundance gets larger 
as we decrease the spin-independent elastic scattering cross 
section leading to the diagonal edge clearly visible in the dis- 
tribution of points. 

The second region corresponding to DM co-annihilations 
in FIG. |4]) has generically lower spin-independent scatter- 
ing cross sections but can have a range of relic abun- 
dances. The majority of points in this region correspond 
to situations where the freeze-out process is unrelated to 
the spin-independent cross section as is the case for stau 
co-annihilations and exchange and so no discernible pattern 
emerges. 

In the left hand panel of FIG.[4]the relic abundance is kept at 
the freeze-out value and the resulting relaxation of the elastic 
scattering cross section bound is once again apparent due to 
the reduction of the DM relic density compared to the WMAP 
observed value. In this scenario the limits from dSphs also 
play no role due to the suppression in the DM relic density. 

In the right hand panel of FIG. |4j with the DM density re- 
generated to the WMAP observed value, a significant number 
of points are excluded by direct detection. The effect of the 
limits from dSphs is quite minimal, only a handful of points 
(red points in FIG. |4]) are ruled out exclusively by this indi- 
rect constraint and they are the ones with very low freeze-out 
relic abundance and hence a large DM annihilation cross sec- 
tion. In particular, these points represent on-shell annihilation 
through A^. 

FIG.|4]is particularly interesting as it shows that unless the 
cross section is very suppressecj^(asi <C 10~^^cm^), neutrali- 
nos with a freeze-out relic density that exceeds one percent 
of the WMAP upper limit are the only possible type of DM 
candidates that can be saved via a regeneration mechanism. 

In FIG. [5] the same points are shown on plots in the 
(<I>PP,asi) plane. These plots give a useful demonstration of 
the relative importance of the two constraints, with the major- 
ity of points being ruled out by direct detection. 



^ Even with this suppression the number of points in this region is very low. 
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FIG. 6. A plot of m^^o versus the freeze-out neutralino relic den- 
sity where regeneration is assumed. Colour coding is the same as in 
FIG. [3] 

FIG. [6] shows the final result of applying both direct and 
indirect detection constraints in the (m^o, ^Ipo^^) plane as- 
suming the regeneration of the DM relic abundance to the 
WMAP observed value. It can be seen that indirect detection 
limits do not constrain the resonant Z and freeze-out anni- 
hilation scenarios. Spin-independent direct detection excludes 
the most under abundant scenarios particularly in the case of 
resonant annihilation via . The interplay between the spin- 
independent coupling and resonant effects during freeze-out 
discussed earlier can again be seen in the thin strip of points 
excluded around the edges of the Z and resonance regions. 
It is clear that points further from the resonance regions re- 
quire larger couplings in order to reduce the freeze-out relic 
abundance below the WMAP observed value. This generates 
a larger spin-independent cross-section leading to the exclu- 
sion of these points by direct detection. 



B. Regeneration in scenarios with m^o > 100 GeV 

Xi 

Turning now to the heavier candidates of Scan B. FIG. [7] 
displays the distribution of MCMC points found in the 
(<I>PP,m^o) plane. The colour scheme is identical to the earlier 

Xi 

figures with red points ruled out by the dSph limits, yellow 
points ruled out by the XENON 100, grey points ruled out by 
both and green points are not constrained by either. 

The left panel of FIG.|7j with no regeneration, has no points 
that are ruled out by the dSph limits. As with neutralinos with 
masses below 100 GeV, the dSph limits plays no significant 
role in restricting the under- abundant scenarios due to the re- 
duced relic density suppressing the DM annihilation rate into 
photons. 

Regenerating the DM density to the WMAP observed 
value, the dSph limits now play a significant role in constrain- 
ing the allowed parameter space as demonstrated in the right 
panel of FIG. [7] The points ruled out by the dSph limits cor- 
respond to the most under- abundant scenarios, which can be 
seen clearly in FIG. [8j which contains plots of points in the 
(Gsi, ^^FO^^) plane. In the case of regeneration (right panel of 
FIG. [8]) the impact of the dSph limits is restricted to the most 
under-abundant scenarios with abundances up to just below 
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FIG. 7. <I>pp vs the neutralino mass for Scan B. The Hmit shown is from FERMI-LAT observations of dSphs ^9\. The right panel shows the 
case with the regeneration of the DM rehc density to the correct value, the left panel shows the case without. Colour coding is the same as in 
FIG.|3l 
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FIG. 8. Spin-independent cross section versus the neutralino relic density with m^o > 100 GeV. The right panel shows the case with the 
regeneration of the DM relic density to the correct value, the left panel shows the case without. Colour coding is the same as in FIG. [3] 




m-o [GeV] 



FIG. 9. Neutralino mass versus the freeze-out neutralino relic density 
where the regeneration of the DM density is assumed. Colour coding 
is the same as in FIG. [3] 



3% of the WMAP observed value being constrained. 

In addition, FIG. [8] shows that direct detection still plays 
an important role in constraining neutralino DM with masses 
above 100 GeV. In particular, it constrains points with a large 
range of freeze-out abundances and consequently provides a 
useful complementary constraint to the dSph limits. 

The result of applying the indirect and direct detection lim- 
its in the (Q^^oh^ ,m^o) plane after regeneration is shown in 
FIG. [9] The majority of excluded points come from the lower 
end of the mass distribution with all points with relic abun- 



dances less than around 3% of the WMAP observed value 
being ruled out by a combination of the direct detection and 
dSph limits applied in our analysis. Spin-independent direct 
detection limits also lead to a reduction of points with larger 
masses and abundances. 



IV. DISCUSION AND CONCLUSION 

Using a familiar framework (neutralinos in the MSSM), we 
have investigated the configurations for which the expected 
freeze-out relic density could be much smaller than the upper 
limit of the WMAP observed value. We have found many con- 
figurations where ^Ipo^^ could be down to 10~^^1wmap^^- 
In particular, for low neutralino masses, resonant annihilation 
through Higgs or Z boson appear to be very efficient. 

However, we have demonstrated that if a mechanism is ca- 
pable of regenerating the candidate DM number density to 
the present observed value, the combination of FERMI-LAT 
gamma ray observations in dSph and DM direct detection lim- 
its from XENON 100 make such scenarios difficult to realise, 
thereby suggesting that candidates with very small freeze-out 
relic density (less than a percent of the WMAP upper limit) 
cannot be the sole explanation to the DM problem even if one 
assumes that after freeze-out the DM density is regenerated. 

One of the central points of this study is the link between 
the DM annihilation process at freeze-out and the predicted 
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direct detection rate. Essentially, as is already known in the 
MSSM scenario a light Higgs or Z exchange is needed in or- 
der for light neutralinos to become extremely under-abundant. 
Small neutralino freeze-out abundances correspond to scenar- 
ios which are close to the Higgs resonance and thus typi- 
cally predict large neutralino-nucleon elastic scattering cross 
sections since one is close to the minimal Higgs mass value 
that is possible in the MSSM (given collider constraints) 
and the neutralino-Higgs couplings are constrained by the re- 
quirement of a large annihilation cross section which is itself 
bounded by the Higgs decay widtlj^ 

In scenarios without regeneration, such a large neutralino- 
nucleon elastic scattering cross section (or large annihilation 
cross section) is not necessarily excluded. Indeed, the very 
small freeze-out abundance actually induces a 'r|' suppression 
factor and reduces the elastic scattering rate in nuclear recoil 
direct detection experiments as well as the indirect detection 
rate. However in scenarios where one allows regeneration of 
the relic density to happen, there is no 'r|' suppression factor 
and these scenarios can be ruled out by direct and indirect 
detection experiments. 

The heavy neutralino scenarios which are under- abundant 
(less than 3% of the observed relic density) also benefit from 
the 'r|' factor suppression if there is no regeneration mecha- 



nism involved, so they cannot be ruled out. However, when 
we assume regeneration, we find that the indirect detection 
constraint set by the FERMI-LAT experiment actually rule 
out these candidates and complement the constraint set by 
XENONIOO. 

Our conclusion is based on the combination of astrophysics, 
astroparticle and particle physics data. Any more constraints 
in, at least, one of these fields will enable stronger constraints 
to be set, thus restricting the types of mechanisms that could 
give the DM the relic density it has today. 

In this analysis we have not applied the latest limits emerg- 
ing form the LHC on the sparticle spectrum. Applying these 
limits will reduce the number of allowed points further mak- 
ing the regeneration scenario even harder to realise. We leave 
this study to future work. 
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